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Abstract 

The phase oscillator model with global coupling is extended to the case of finite- 
range nonlocal coupling. Under suitable conditions, peculiar patterns emerge in 
which a quasi-continuous array of identical oscillators separates sharply into two 
domains, one composed of mutually synchronized oscillators with unique frequency 
and the other composed of desynchronized oscillators with distributed frequencies. 
We apply a theory similar to the one which successfully explained the onset of 
collective synchronization in globally coupled phase oscillators with frequency dis- 
tribution. A space-dependent order parameter is thus introduced, and an exact 
functional self-consistency equation is derived for this quantity. Its numerical solu- 
tion is confirmed to reproduce the simulation results accurately. 
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1 Introduction 



Large populations and continuous fields of coupled oscillators form a represen- 
tative class of synergetic systems met in a wide range of scientific disciplines 
from physics, chemistry, engineering, biology to brain science [|lpf| • Collective 
dynamics of coupled oscillators depends crucially on the range of their mutual 
coupling. It has recently been realized that when the coupling is nonlocal, 
the patterns which emerge could be drastically different from those which we 
expect for oscillators with local or global coupling The implication of 

this fact is relevant even to what we conventionally call locally coupled sys- 
tems, typically reaction-diffusion systems. This is because it may happen that 
nonlocality can arise effectively as a result of elimination of some variables, 
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e.g., rapidly diffusing components in the case of reaction-diffusion dynamics. 
Among the variety of patterns which are characteristic to nonlocally coupled 
oscillators, we will focus our attention below on a particular class of patterns 
in which the whole medium is separated into two domains of qualitatively 
different dynamics. Specifically, the oscillators are mutually synchronized in 
one domain while they are completely desynchronized. in the other domain. A 
preliminary work on such dynamics was reported recent lyJ8|. We will present 
below a more thorough investigation of this problem. 

The collective dynamics of our concern is similar to the collective synchro- 
nization in globally coupled oscillators with distributed natural frequencies 
PJ9|JT0[| where the whole population splits into two subpopulations each com- 
posed of synchronized and desynchronized oscillators. There the systems is 
stationary in a statistical sense within a constant drift of the collective phase 
corresponding to the oscillation of the population as a whole. Similarly to the 
case of global coupling, the individual oscillators can be regarded as being 
controlled by a common mean field, although the latter is now space depen- 
dent. The oscillators under consideration are identical in nature, still we may 
develop an exact mean-field theory similar to the one developed for globally 
coupled oscillators. Commonly to these theories, a suitably defined order pa- 
rameter representing the collective state is determined from a condition for 
self-consistency which must exist between the collective dynamics and the 
dynamics of the individual oscillators. 

In section 2, some results of our numerical simulation will be presented. The 
simulation is carried out first on a nonlocal version of the complex Ginzburg- 
Landau (CGL) equation in one space dimension with periodic boundary con- 
ditions, and then on another equation obtained through the phase reduction of 
the nonlocal CGL. In section 3, our theory will be developed and a functional 
self-consistency equation will be derived for the order parameter. Then the 
numerical solution of the self-consistency equation will be compared with the 
simulation results presented in Section 2. Concluding remarks will be given in 
the final section. 



2 Coexistence of coherence and incoherence 

As a simplest model for densely and uniformly distributed oscillators with non- 
local coupling, let us consider the following equation for a complex amplitude 
A which we call nonlocally coupled complex Ginzburg-Landau equation^. 
In one spatial dimension our model is expressed as 

d 

TrA(x,t) = (l + iu )A- (l+ib)\A\ 2 A + K(l + ia)(z(x,t) - A{x,t)).{l) 
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The quantity Z(x,t), which we call the mean field, represents the effects of 
the nonlocal coupling and is given by 



Z(x, t) = J G(x - x')A(x', t)dx. (2) 
The coupling function G changes with the distance as 

G (y) = 2 ex P( -K M) ( 3 ) 

and is normalized. The exponential form of G is in fact a natural consequence 
of the reductive derivation of Eq. (1) from a certain class of reaction-diffusion 
systems when the latter involve an inactive diffusive component to be elimi- 
nated adiabatically||[7|]. Note that the last term in Eq. (1) representing the 
coupling is so arranged that it may vanish when A is uniform in space. If A is 
nonuniform but sufficiently long-waved, this term can be approximated with 
a diffusion term, and then Eq. (1) is reduced to the ordinary CGL. 

Equation (1) is numerically simulated on a finite interval [0, 1] with periodic 
boundary conditions. In a suitable range of parameters, in which the uniform 
oscillation is linearly stable, the system can develop peculiar patterns when 
the initial perturbation is finite. Figure 1 shows a typical example of such 
phase patterns. The initial condition is such that the modulus of A equals 1 
everywhere and its phase changing randomly still its envelope being nearly 
symmetric about the midpoint x = 1/2. Since the coupling constant K is as- 
sumed to be relatively small, the local oscillations attain almost full amplitude, 
so that we will mainly be interested in the phase pattern. The pattern consists 
of two distinct domains. In one domain which appears near the boundaries, 
the pattern is spatially continuous and smooth, whereas in the central domain 
spatial continuity seems completely lost. The pattern as a whole looks steadily 
advancing upward, implying that the entire population is undergoing regular 
collective oscillation with a definite frequency Q. 

To investigate the nature and origin of such patterns in further detail, Eq. 
(1) is reduced to a phase equation which would be much easier to analyze. 
The phase reduction is valid when the coupling strength K is small, which is 
actually the case, and the reduced equation takes the form 

t) = uj — J G(x — x') sin(0(x, t) — 4>{x\ t) + ajdx', (4) 

where the time scale has been so changed as to normalize the coupling strength, 
u is the natural frequency rescaled accordingly, and the phase constant a in 
the phase-coupling function is related to the original parameters through 

b — a , ,_ x 

tana = -, a(b — a) > 0. (5) 

1 + ah y ' w 
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Under the parameter condition corresponding to the previous numerical sim- 
ulation of Eq. (1), a phase pattern similar to FIG. 1 was obtained, which is 
displayed in FIG. 2. 

We now focus on the phase patterns obtained for our phase equation. Let us 
introduce relative phase ip by 

= at + ip (6) 

which describes the dynamics of the phase deviation from the reference motion 
with some drift velocity Q whose value is still open. We rewrite Eq. (4) using 

ip as 

^(x, t)=u-Sl- J dx'G(x - x') sin(V>(z, t) - ip(x', t) + a) . (7) 

As a generalization of the theory of synchronization transition of globally 
coupled phase oscillators with frequency distribution, we introduce a complex 
order parameter with modulus R and phase through 

J dx'G(x - x') exp[i^(x' ,t)} = R(x,t) exp[i&(x,t)]. (8) 

Unlike the order parameter defined for globally coupled oscillators, the above 
quantity is space-dependent. Still a physical picture is valid such that we 
are practically working with an assembly of independent oscillators under the 
control of a common forcing field represented by R and 0. This is confirmed by 
rewriting Eq. (7) in terms of the order parameter into a forced one-oscillator 
equation as 

d 

— ip( x , t) =cu-tt- R(x, t) sm(^j(x, t) + a + 0(x, tfj . (9) 

The spatial profiles of R(x) and Q(x) after a long-time average are displayed in 
FIG. 3a and 3b, respectively. We see that the forcing mean-field amplitude is 
stronger near the boundaries and weaker near the center of the system. If the 
phase distribution such as shown in FIG. 2 is statistically stationary within 
a constant overall drift with velocity Q, then R and are time- independent. 
Their approximate constancy is confirmed by the fact that the fluctuation of 
the time sequence of R(x,t) observed at the midpoint of the system x — 1/2 
gives a standard deviation as small as 0.008 which comes possibly from a 
finite-size effects. Similar property can also be confirmed for Q(x,t). 

Figure 3c shows a distribution of the actual frequencies cu(x) of the individual 
oscillators. This was also obtained through a long-time average. In the coherent 
domain, the oscillation frequencies have an identical value Q, while in the 
incoherent domain they are distributed but give a well-defined continuous 
curve. 
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The dynamics which is going on is now clear. The system is divided into 
two subgroups of oscillators. In the first group, the forcing amplitude is large 
enough for the oscillators to be entrained, so that they oscillate with an iden- 
tical frequency Q. In the second group, in contrast, the forcing amplitude is 
too weak for entrainment, so that the frequencies of the individual oscillators 
differ from Q. Since the actual frequencies of these desynchronized oscilla- 
tors should depend on the local amplitude of the forcing, the nonuniform but 
smooth pattern of R(x) as shown in FIG. 3a implies a smooth distribution 
of the actual frequencies, and this is consistent with FIG. 3c showing u(x). 
Our next problem is how to predict stationary patterns of R(x) and Q(x) 
theoretically, and this will be the subject of the next section. 



3 Theory 

The basic equations to work with below are Eqs. (8) and (9). Assuming R and 
to be time-independent, we first try to find solutions of Eq. (9) for the phases 
ip(x,t) each of which should be a function of R(x) and Q(x). These solutions 
are then substituted into Eq. (8), leading to an equation for determining R and 
as a functional of R and in a self-consistent manner. Since the collective 
frequency Q is still to be determined, we will be working with a nonlinear 
eigenvalue problem. 

Here a question may arise as to the seeming contradiction between the as- 
sumed time-independence of R and and the apparent time-dependence of 
the solution of Eq. (9) which is true if R(x) is below a certain critical value. 
The same problem was encountered in the problem of synchronization tran- 
sition in globally coupled phase oscillators with frequency distribution. The 
answer to this question is that the factor exp[iip(x', t)] appearing in the inte- 
gral of Eq. (8) should be replaced with is statistical average (exp[iif)(x' ,t)]). 
Such an average can be calculated using a suitable invariant measure which is 
easy to find. 

The successive steps to be taken for reaching the final solution are the follow- 
ing. Note first that Eq. (9) admits two types of solutions depending on the 
location of the oscillators of concern. They are constant solutions and drifting 
solutions. If \R/{oj — fi)| < 1, we have a time-independent solution given by 



If \R/{oj — O) | > 1, in contrast, we have a drifting solution. The actual fre- 
quency cu of the corresponding desynchronized oscillator as a function or R is 




(10) 
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given by 



uj = Q + (cu — Q) > 



which equals the collective frequency Q plus the mean drift velocity of ip. 
Instead of using such drifting solutions in Eq. (9), we use an invariant measure, 
i.e., the probability density p(9, R) such that ip takes on value 9 for given R, 
which is inversely proportional to the drift velocity of ip at ip — 6. Thus, 

p(6,R) = C(u-n-Rsm6)- 1 , (12) 



where C is a normalization constant given by C = (2ir)~ y (u — Vt) 2 — R 2 . 
As noted earlier, this invariant measure is used for replacing the exponential 
factor appearing on the left-hand side of Eq. (8) with its statistical average. 
The physical idea underlying this replacement is that an equilibrium statistical 
ensemble may be applicable to each small local subsystem because infinitely 
many oscillators which are independently oscillating are already contained in 
such a subsystem. 

The spatial profile of R(x) in FIG. 3 implies that there are spatial points 
x = \ ± x c at which R takes a critical value separating the two types of 
solutions described above. Namely, in the inner domain (\x' — || > x c ) the 
solutions are drifting, while in the outer domain (\x' — || > x c ) they are 
constant in time. Thus Eq. (8) becomes 



i 



J dx'G(\x-x'\)ex.p[i(0(x') - a)]h(R(x'j) = R(x)exp[iO(x)], (13) 
o 

where 

J* n exp(i9)p(9,R)d9 (\x'-±\< x c ) 

The above expression for h(R) in each domain can further be made explicit 
as 



expK(«)] = /^^) I + ^, (15) 

J exp{i0)p(6, R)d6 = ^[uj-VL- ^(uj- Q) 2 - R 2 }. (16) 
o 

Equations (13)~(16) constitute a functional self-consistency condition from 
which R(x), Q(x) and Q are to be determined. By using R(x) and fl thus ob- 
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tained, the actual frequencies of the desynchronized oscillators are determined 
from Eq. (11). 

There is a unique value Q of the collective frequency for which our functional 
self-consistency equation admits a solution. Such solution can be obtained 
numerically, and is compared with our numerical simulation. Some of the 
results are given in FIG. 4a~4c. We find that the agreement between the 
theory and numerical experiments is almost perfect for each quantity of R(x), 
Q(x), Q and cu(x). 



4 concluding remarks 

The present theory makes full use of the fact that the number of oscillators 
contained in the range of coupling is practically infinite. This fact allows us 
to apply a mean field theory, and our problem is reduced to a one-oscillator 
problem supplemented with a self-consistency condition. It is well known that 
the same idea works nicely for globally coupled oscillators. The present theory 
suggests that, beyond a special problem discussed in the present paper, the 
mean-field idea may work more generally for nonlocally coupled systems as far 
as the coupling radius remains much longer than the minimum length scale 
associated with the discreteness of the spatial distribution of the constituent 
oscillators. 
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Figure Captions 



FIG.l. Instantaneous spatial distribution of the phases obtained from Eq. (1) ap- 
proximated with an array of 512 oscillators distributed over an interval of 
length 1 with periodic boundary conditions. Parameter values are a = —1.0, 
b = 0.88, k = 4.0. 

FIG. 2. Instantaneous spatial distribution of the phases obtained from Eq. (4) ap- 
proximated with an array of 512 oscillators distributed over an interval of 
length 1 with periodic boundary conditions. Parameter values are a = 1.457, 
k = 4.0. 

FIG. 3. Spatial profiles of the long-time average of the order parameter amplitude 

R(x) (a), the order parameter phase Q(x) (b) and actual frequencies Q(x) 

(c), each obtained from Eq. (4). 
FIG. 4. Theoretical curves for R(x) (solid line in (a)), Q(x) (solid line in (b)) and 

Q(x) (broken line in (c)) are compared with the corresponding quantities in 

FIG. 3 obtained numerically. 
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